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Problem statement 


unique Tb observation —Xunique SM retrieval 
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Problem statement 


unique Tb observation —Xunique SM retrieval 


Because: 
¢ choice in RTM parameterization (e.g. SMAP L2, LMEB L2, ECMWF, SMAP L4) 
¢ choice in inversion algorithms: 

o regular (i.e. non-mpdi-based) or mpdi-based algorithm 


o species included in cost function (CF): H-pol, V-pol, which angle(s), how 
many angles? 


o whether to include proir soil moisture information in the CF 
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Problem statement 


unique Tb observation —Xunique SM retrieval 


Because: 
¢ choice in RTM parameterization (e.g. SMAP L2, LMEB L2, ECMWF, SMAP L4) 
¢ choice in inversion algorithms: 

o regular (i.e. non-mpdi-based) or mpdi-based algorithm 


o species included in cost function (CF): H-pol, V-pol, which angle(s), how 
many angles? 


o whether to include proir soil moisture information in the CF 


> note: uncertainty = systemic error + random error 
random error may be focus (e.g. in data assimilation studies) 


KU LEUVEN 


J. Quets et al. (KULeuven) Uncertainty Assessment of Space-Borne Passive Soil Moisture Retrievals SMW 2017, Vienna 


11 EASEv2 grid cells containing SMAP core validation sites 


Goole Earth 
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11 EASEv2 grid cells containing SMAP core validation sites 


Little River watershed 
36 km EASEv2 grid 
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> in situ soil moisture observations to which SM retrievals will 
be compared (May 2010 June2015) 
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Methods: ensemble sets 


¢ 4choices in RTM parameterization 


o Lit1: based on SMAP L2 h: 0.1140.16; w: 0.0530.07; b,&b,: 0.130.11 
o Lit2: based on LMEB L2 h: 0.100.70; w: 0.050.05; b,&b,: 0.1540.3 
o Lit3: based on ECMWF h: 1.661.66; w: 0.0030.05; b,&b,: 0.1530.3 
o Lit4: based on SMAP L4 h: 0.000.97; w: 0.0030.13; b,&b,: 0.0740.4 


¢ 4perturbations for each h,,,, Nya, W, and b,, b, (-50%, -25%, +25%, +50%) 

¢ 7 angles in CF (i.e. 30°, 35°, 40°, 45°, 50°, 55°, 60°), either separately or together 
* 2 polarizations (i.e. H-pol, V-pol) 

¢ 2 different RTM-inversion algorithms (i.e. mpdi-based or non-mpdi-based) 
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Methods: ensemble sets 


¢ 4choices in RTM parameterization 


o Lit1: based on SMAP L2 h: 0.1140.16; w: 0.0530.07; b,&b,: 0.1 30.11 
o Lit2: based on LMEB L2 h: 0.1030.70; w: 0.0540.05; b,&b,: 0.1530.3 
o Lit3: based on ECMWF h: 1.66 1.66; w: 0.0030.05; b,&b,: 0.1530.3 
o Lit4: based on SMAP L4 h: 0.000.97; w: 0.0030.13; b,&b,: 0.0730.4 


¢ 4 perturbations for each h,,,, Nya, W, and b,, b, (-50%, -25%, +25%, +50%) 

¢ 7 angles in CF (i.e. 30°, 35°, 40°, 45°, 50°, 55°, 60°), either separately or together 
¢ 2 polarizations (i.e. H-pol, V-pol) 

¢ 2 different RTM-inversion algorithms (i.e. mpdi-based or non-mpdi-based) 

* +> many ensemble sets tested, ranging between 28 and 2456 members 

¢ +> not all combinations possible 
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Methods: Cost Function (CF) 


=s 1 
CF = (TDsim — TDops)'C . (TDsim TDops) si 0,022 (SM -etr _ SM cism)* 


with C = Tb error covariance matrix, representing: 
- Tb error variances (62 K?) 
- correlations between Tb errors of different incidence angles 


with prior SM information included 
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Results: retrieval uncertainty 


Part 1: sensitivity analysis 
retrieved soil moisture 


default retrieval: - 7 site: Forth Cobb 
| ae =o 
° single species in CF: 40° Hpol = 0.2 | 
‘ . . . — l i IN al NA ! LAW I ! Als 1 
Lité RTM parameterization = "3010 2011 2012 2013 2014 2015 
*  non-mpdi-based inversion algorithm ” time 


° =basically SCA 
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Results: retrieval uncertainty 


Part 1: sensitivity analysis 
retrieved soil moisture 


default retrieval: site: Forth Cobb 
al 
* single species in CF: 40° Hpol 0.2 | <a 
é . . . — L lt IN al NA ! LAW UM i i Als 1 
Lit4 RTM parameterization = *010 2011 2012 2013 2014 2015 
° non-mpdi-based inversion algorithm “” time 


° =basically SCA 
HOW? 


° step1: choosing angle, polarization, RTM-parameters, inversion algorithms separately 
° step2: calculating ensemble variances of these experiments 


° step3: dividing this variance in long-term mean ensemble variance and short-term ensemble variance 
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Results: retrieval uncertainty 


Part 1: sensitivity analysis 
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Results: retrieval uncertainty 


Part 1: sensitivity analysis 
retrieved soil moisture 
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Results: retrieval uncertainty 


Part 1: sensitivity analysis 
retrieved soil moisture 
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Results: retrieval uncertainty 


Part 1: sensitivity analysis 
retrieved soil moisture 
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retrieved soil moisture 
site: Forth Cobb 


inversion algorithm: regular 


inversion algorithm: mpdi-based 
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Results: retrieval uncertainty 


Part 2: total uncertainty estimation 
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Results: retrieval uncertainty 


Part 2: total uncertainty estimation 
> (1) find a properly verified ensemble set 
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Results: retrieval uncertainty 


Part 2: total uncertainty estimation 
> (1) find a properly verified ensemble set 
> (2) its ensemble variance characterizes total retrieval uncertainty 
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Results: retrieval uncertainty 
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Results: retrieval uncertainty 
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Results: retrieval uncertainty 
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Results: retrieval uncertainty 
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Results: retrieval uncertainty 


Part 2: total uncertainty estimation ensemble verification 
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Results: retrieval uncertainty 


Part 2: total uncertainty estimation ensemble verification 
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time-averaged ensemble variance of verified ensemble set =78% 
of variance of jn situ observations 
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Results: ranked skills of ensemble retrievals 
Ranked Skills (n=2856) 
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Results: skills of ensemble means 


Skill metrics of ensemble means VS ensemble sets 
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Results: skills of ensemble means 


Skill metrics of ensemble means VS ensemble sets 
HE ensemble mean 

operational SMOSV6 

CLSM 


no CLSM constraint in CF CLSM constraint in CF 


0.06 
ase | 
2 


bias: 


ubRMSE: 


bias: 


> inclusion of CLSM 
sm in the CF generally 
improves every skill 
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Results: skills of ensemble means 


Skill metrics of ensemble means VS ensemble sets 
HE ensemble mean 
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> inclusion of CLSM 
sm in the CF generally 
improves every skill 

> especially when 
including all species in 
the CF 
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> inclusion of CLSM 
sm in the CF generally 
improves every skill 

> especially when 
including all species in 
the CF 
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Results: skills of ensemble means 


Skill metrics of ensemble means VS ensemble sets 
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Take home messages 


1. passive L-band SMOS soil moisture retrievals are uncertain 


... and most sensitive to RTM parameterizations (e.g. roughness parameters and 
surface albedo) 


... with the ensemble variance of a verified set amounting to 78% of in situ temporal 
variance 


... choice of RTM-parameter set strongly influences the bias 

2. constraining a CF with CLSM-simulated soil moisture improves the retrieval skill 
even though CLSM skills are generally worse than retrieval skills 
main reason: constrain extreme high and low values 

3. ensemble means of ensemble sets 


ensemble means of ensemble sets outperform operational SMOS by about up to 
9% for UDRMSE and more than 6% for anomaly R 


best performance reached by including as many as possible species in the CF (i.e. 
14 species) 


4. next: compare to SMOS-IC or other alternatives 
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